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1. Introduction 

1.1. Hermitean Wishart random matrices: A brief overview 

Empirical covariance matrices C is a central object in statistical analysis of dynamical 
complex systems characterised by multivariate time series. Built upon n discretised 
sets of 'experimental' signals {xi(t), . . . , x n (t)} of the length m each[f] 

m 

<w = y^ x a(tj) Xgi(tj), (a,a') G (1, ...,n), (1.1) 

3=1 

these matrices are intrinsically noisy. The experimental data (i a (t)} can be as diverse 
as (i) time-dependent price changes of various assets (Laloux, Cizeau, Bouchaud and 
Potters 1999; Plerou, Gopikshnan, Rosenow, Amaral and Stanley 1999; Plerou, Gopik- 
shnan, Rosenow, Amaral, Guhr and Stanley 2002) at a stock exchange; (ii) multichan- 
nel physiological (e.g., electro- and magnetoencephalography) recordings (Kwapich, 
Drozdz, Liu and Ioannides 1998; Kwapieh, Drozdz and Ioannides 2000; Seba 2003); 

(iii) an atmospheric parameter (such as wind velocity, geopotential height, or temper- 
ature) as a function of time and a discretised coordinate (Santhanam and Patra 2001); 

(iv) information flows in the world-wide web (Barthelemy, Gondran and Guichard 
2002); (v) pixel decomposed signal of a noisy image (Basu, Ray and Panigrahi 2010), 
and many more. 

The random content of empirical covariance matrices originates from (i) finiteness of 
the time series used (when the length m of a time series is not very large compared 
to n, the number of experimental signals) and (ii) non-stationarity of system corre- 
lations (Plerou, Gopikshnan, Rosenow, Amaral and Stanley 2000). As this intrinsic 
noise blurs the information about genuine correlations in a complex system, empirical 
covariance matrices must be 'cleaned'. A pretty efficient denoising can be achieved by 
confronting statistics of eigenvalues and eigenvectors of empirical covariance matrices 
C with those of most random covariance matrices W constructed from mutually uncor- 
rected (Gaussian) times series. Evidently, the deviations in statistics of C and of W 
are to feature the properties that are specific to a complex system and can therefore be 
used to extract genuine correlations between various system components out of experi- 
mental data x a (t). In this paper, we shall only be concerned with the spectral analysis 
of covariance matrices. For an introductory exposition of the eigenvector-based analy- 
sis the reader is referred to the review by Bouchaud and Potters (2009). 

The most random covariance matrices, also known as the Wishart matrices (Wishart 
1928), are fundamental to the whole realm of multivariate statistical analysis (Muirhead 
1982). Given n sets of uncorrelated discretised zero- mean Gaussian random processes 
X a (i), an n x n random Wishart matrix W is defined by its entries 

m 

W aa > =Y2 X a(tj)X<*'(t j ), (a,a r ) G (l,...,n). (1.2) 

3=1 

In a matrix notation, 

W = XX\ (1.3) 

where X denotes a rectangular nx m matrix whose (a,j) entry equals X a .j = X a (tj), 
and m is the length of a time series. The associated probability measure reads 

dPi%(X) = c^ m (P,a ) exp [-aptrXX^] DX. (1.4) 

X The a-th signal x a (t) is a vector consisting of m components: x a (i) = {x a (ti), . . . ,x a (t m )}. 



1 Introduction 



4 



Here, c n , m {P) is the normalisation constant, 

v Bnm/2 
7T \ 



Cn, m (P,CLp) = , (1.5) 



the notation Z?Af stands for 

n m f fi 

As) 



^=nn n< ■ 

a=l j=l \s=l / 



and /? denotes the Dyson index (Mehta 2004). The superscript (s) in Eq. (1.6) counts 
a number of degrees of freedom in X a (tj). For real, complex and quaternion real time 
series, /3 equals 1, 2, and 4, respectively. Equations (1.3 1 - (1.6 1 uniquely define three 
canonical ensembles of random Wishart matrices. 

By definition, the canonical Wishart probability measure dir^m(W) is provided by 
the matrix integral 

d*%L0V) = J dPi%(X) 6 (w - ATAft) DW ; (1.7) 

where the flat measure DW equals 

n n / P \ 

a— 1 a>a' — 1 \s — 1 / 

In the Wishart domain (to > n), this matrix integral can be evaluated in a closed form 
resulting in the elegant formula (Muirhead 1982) [|] 

dnglm = c-l(l3,ap)G(W) (det W) /3(n -" l+1)/2 " 1 exp[-a^trW] DW. (1.9) 
Here, c„ iirl (/3, a^) is yet another normalisation constant, 

n(n-l)/3/4 " / fl \ 

<*) = nnW2 n r ?( m + j - !) • 

The Heaviside function 6(W) indicates that W is a positive definite matrix. 

Spectral statistical properties of the Wishart random matrices are well understood 
since the probability measure in the eigenvalue space readily follows from Eq. (1.9| 
upon diagonalisation of W (Mehta 2004) : 

n n 

dp^(Xi, . . . , A n ) = c-i m (P, a,) J] e(A,) A J(™-«+D/2- |A n (A)|^ J] d\j. 

(1.11) 

Here, c n ,m(P,a>p) is the normalisation constant, 

«.) = ar /. [r ; i + , /2) i,. fir (i + ^) r (i + f (m - » + i)) . (1.12) 

Spectral correlation functions of arbitrary finite order can be restored from Eq. ( |1.11| 
by applying the orthogonal polynomial technique (Mehta 2004). The mean density 
of eigenlevels Q^) m {X) is of primary importance for denoising of empirical covariance 

§ The anti- Wishart domain (m < n) has been treated by Janik and Nowak (2003). 
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matrices. Its large- n limit (when the ratio q = m/n > 1 is kept fixed) is described by 
the Marcenko-Pastur law (Marcenko and Pastur 1967) 



A) (A — A m i n ), 



where A € (A min , A ma > 



f3n A 

). The endpoints of the eigenvalue support are 



A {.™}^2^ ( ^ Tl) 



(1.13) 



(1.14) 



Equations (1.13) and (1.14) lay the basis for a cleaning procedure required to separate 



information-carrying correlations from the useless noise which is always present in 
empirical covariance matrices (Bouchaud and Potters 2009). 

1.2. Non-Hermitean Wishart random matrices 
1.2.1. Motivation and definitions 

Empirical covariance matrices C discussed in the previous Section have a built-in Hcr- 



miticity [Eq. (1.1)]. However, in many applied problems of statistical analysis the 
empirical covariance matrix can equally be non-Hermitean. This is the case when one 
is interested in studying correlations between two 'remote' complex systems charac- 
terised by two distinct sets of time series {aji(t), . . . , x n (t)} and {y 1 (t), . . . , y n (t)}, re- 
spectively. These series may represent multichannel magneto-encephalography record- 
ings of activity in the left and right auditory cortex (Kwapien, Drozdz and Ioannides 
2000) or describe return of stocks traded on two large but geographically distant mar- 
kets (Kwapien, Drozdz, Gorski and Oswiccimka 2006). In both cases, a noise-dressed 
empirical covariance matrix 



r , — 



3=1 



x a (tj)y a '(tj), (a,a r ) e (l,...,rc), 



(1.15) 



is manifestly non-Hermitean. Consequently, its spectrum as well as the spectrum of an 
appropriate most random covariance matrix 



W aa > 



j=l 



X a (tj ) Y a i [t 



j )• 



(a, a') e (1, 



(1.16) 



composed of 2n sets of uncorrelated discretised zero-mean Gaussian random processes 
X a (t) = {X a (ti),...,X a (t m )} and Y a (t) = {Y a (tx), . . . , Y a (t m )}, becomes com- 
plex valued. The latter observation has been beautifully illustrated in the studies by 
Kwapien, Drozdz and Ioannides (2000) and Kwapien, Drozdz, Gorski and Oswiccimka 
(2006). 

This reasoning leads us to introduce three non-Hermitean counterparts of canonical 

(1.17) 



Wishart matrix models [Eqs. ( |1.3| and (1.4)] defined as 
The probability measures assigned to the matrices X and y read 



dP(%(X) = c~) m (P,ap)e X p 



-a/3 tr XX 



dPi%(y) = c-UMb) exp \-aUryy 



DX, 

Dy. 



(1.18) 
(1.19) 



The normalisation constant c„ im (/3,a^) and the flat measures DX and Dy are speci- 
fied by Eqs. (1.5) and (1.6), respectively. 



The Dyson fluid approach is presumably the shortest way to derive this result, see Dyson (1971). 
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Given the above definition of non-Hermitean Wishart random matrices, we would like 
to determine (i) the probability measure induced on W, (ii) the joint probability den- 
sity function of all n complex eigenvalues of W, (hi) the eigenvalue correlation functions 
of arbitrary (finite) order, and also analyse (iv) the mean spectral density in various 
scaling limits. Even though the real valued time series (/3 = 1) are of most interest 
for statistical applications, the present paper will focus on a somewhat simpler case of 
complex valued time series (/3 = 2) [^j] 



1.2.2. Main results and discussion 



For the benefit of the readers, we collect our main results into this easy to read sub- 
section with pointers to the sections containing detailed derivation of each statement. 

(i) The probability measure in the matrix space. — Let W be a matrix drawn from 
an en semb le of n x n co mplex non-Hermitean Wishart random matrices as defined by 
Eqs. (1.171, (1.181 and (1.19). Then, for m > n, the probability measure d7r njm (W) 



associated with the matrix entries {W a . a '} equals 

1 



(Type I) dn n , m (W) = 



(2tt) 2 " 
xdet" ? 



Dq exp 



- tr (qW 



In 



1 



4<22 al. 



where the flat measure DW is defined by 

n 

DVV= Yl dSfitWca' dOmW a , a > 



a,a' — l 



In Eq. (1.20), that will be referred to as the type I representation of d7r r , 



(1.20) 

(1.21) 
,(VV), the 



integration runs over an n x n generic complex valued matrix. 



An alternative, type II representation, 



(Type II) d7r„, m (W) = 



M) n(n + m)/2 n£i n r(j) 



DQQ(Q) 



7r n(3n-i )/ 2 n - =l r(j) y Qt= 

xdet m - 2 "Q exp \-^/a^af 2 tx (Q + VV^Q^W) 



DW. 



(1.22) 

involves an integral over anuxn complex Hermitean matrix. Both type I and type II 
matrix integrals are derived in Section 2.1 and Appendix [A"| 

(ii) The joint probability density function of all complex eigenvalues. — Let (wi, . . . , w n ) 
be complex eigenvalues of an n x n complex non-Hermitean Wishart random matrix. 
Their joint probability density function is 

. 2'> 2 a 2 )"(™+ 1 )/ 2 nr=i n r(i) 
Pn,m\Wx, ■ ■ ■ ,VJ n ) = — 



n" =1 r(i)n7=i r 0') 



ia„mi 2 n 



K 



h\wj\y/a^a/ 2 



(1.23) 



Here, K v {w) is the modified Bessel function of the second kind. For a detailed deriva- 



1 A detailed spectral theory of non-Hermitean Wishart random matrices at /3 = 1 will be a subject 
of a separate publication. 
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Figure 1. Left panel: Numerically simulated distribution of complex eigenvalues 
in the regime {n fixed and v 3> 1}- The red circle of the radius R = is 
displayed to highlight a notable finite-n dispersion effect. Right panel: Profile of 
the mean density Eq. Both figures correspond to n = 10 and v = 95. The 

simulation shown on the left panel involved diagonalisation of 500 matrices. 



tion, the reader is referred to Section |2.2.1| and Appendix [A] Interestingly, the above 
result is a particular case of a more general matrix model introduced by Osborn (2004) 
in the context of quantum chromodynamics (QCD) with a baryon chemical potential 
[see also Akemann, Bloch, Shifrin and Wettig (2008)]. In this connection, let us stress 
that an overlap between our work and the QCD studies cited throughout the paper 
occurs for the nnite-n/fmite-i^ regime. The large-n/large-i^ asymptotic analysis, central 
to statistical applications which triggered our research, is new, and the techniques used 
are different from those employed in the QCD literature. 



(Hi) Determinant structure of spectral correlators. — The p-point correlation function 
follows from Eq. (1.23) by virtue of the Dyson integration theorem (Mehta 1976): 



R nL( w ir ■ ■ ,w P ) = det [K n , m (wj, w fe )] 1 



(1.24) 



where [f] 



fc=0 



2 2k k\(k + v)\ 



is the two-point scalar kernel. Here and everywhere below 

v = m — n > 0. 
See Section 2.2.2 for a straightforward derivation [|] 



(1.25) 
(1.26) 



(iv) The mean density of complex eigenvalues. — An asymptotic behaviour of the mean 
spectral density 

R^(w) =K n<m (w,w) (1.27) 

is of most interest for statistical applications. Depending on the relation between n 
and v (or, equivalently, to, see Eq. (1.26)), the following three limiting laws will be 
established in Section [2~3l 



f To simplify notation, we have set a^a!^ = 1 /4. 

% In view of our previous remark in (ii), Eqs. ( |l.24| and j 1 . 25[ ) are automatically consistent with those 
reported in the QCD literature. 
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Figure 2. Left panel: Numerically simulated distribution of complex eigenvalues 
in the regime (n> 1 and v fixed}. The red circle of the radius R = 2n is displayed 
to indicate the 'mean field' edge of the eigenvalue support (sec discussion in Section 
|2.3,2| , The phenomenon of 'clustering' of eigenvalues around the origin stemming 
from the decay Eq. ( | 1 . 32[ ) is clearly seen. Right panel: Profile of the mean 

density Eq. | |1.29| . A 'goblet base' originates from the erfc-law Eq. | |1.33| . Both 
figures correspond to n = 100 and v = 1. The simulation shown on the left panel 
involved diagonalisation of 100 matrices. 



• Regime I. — For n fixed and v = m, — n 3> 1, we derive the complementary 
T-function law 



This formula assumes that the ratio wj^Jv ~ O(v ) is kept fixed. This result 
mirrors a finitc-n formula for the mean spectral density in the Ginibrc unitary 
ensemble (GinUE) [see, e.g., Eq. (1.43) in Ginibre (1965) and Eq. (2.17) in Ake- 
mann and Kanzieper (2007)]: the two formulae can be reduced to each other upon 
a proper rescaling of the energy variable \w\. 

This regime is characterised by (i) a nearly uniform eigenvalue distribution within 
the disk of the radius R ~ and (ii) a notable (finite-n) dispersion effect that 

manifests itself in a significant number of eigenvalues outside the disk. For an 
illustration, see Fig. [TJ 

• Regime II. — For n 3> 1 and v = m — n > fixed such that v ~ O(n ), a 
complicated behaviour of the mean eigenvalue density is accounted for by the 
single formula 




(1.28) 




1 



I v (\w\)K v (\w\)erfc 



( 



\w\ — 2n 



2v^ 



) 



(1.29) 



It incorporates three different regimes. 



(i) If | it; | ~ O(n ), the above equation reduces to 
R&Hw) = ~I v (\ W \)K u (\w\). 



(1.30) 



For v — 0, the density exhibits a weak logarithmic singularity around the 
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Figure 3. Numerically simulated distribution of complex eigenvalues in the regime 
{n 2> 1 and v 2> 1 with u/n = q kept fixed}. The red circle of the radius 
R = 2n.y / l + q is displayed to indicate the 'mean field' edge of the eigenvalue 
support (see discussion in Section |2.3.3[ l. The phenomenon of 'clustering' of 
eigenvalues around the origin blurs as q increases. Parameters: q = 1/10 (left 
panel), q = 1/2 (middle panel) and q = 1 (right panel). The simulation was 
performed by diagonalising 100 matrices of the size 100 X 100. 



origin; for v > it stays finite: 



2tt \w\ 



={ ' (1.31) 



M<i i I 

— , ^0. 
Air is 



(ii) Increasing 
decay 



to |, one observes that the density Eq. (1.30) exhibits an algebraic 



* u oo 



(w) 



4ir\w\ 



(1.32) 



that holds for \w\ ^> 1. This unusual behavior showing up as a 'clustering' of 
complex eigenvalues around the origin (see Fig. [2]) , contrasts with the uniform 
mean density of complex eigenvalues in GinUE that emerges at \w\ 3> 1 (Gini- 
bre 1965). The 'mean field' result Eq. (1.32) is consistent with a more general 



observation due to Burda, Janik, and Waclaw (2010) who have advocated ex- 
istence of the universal |w| _2 ^ 1_1 / M ^-law for the mean spectral density of the 
product of M independent Gaussian random matrices. We also notice that 
the clustering phenomenon was recently observed in spectra of time-lagged 
correlation matrices (Biely and Thurner 2008). 



(iii) The l/\w\ law breaks down in the y/n- vicinity of the 'critical' point 



2n, 



erfc 



2n 



(1.33) 



where Eq. (1.29) reduces to 

n - ml ' ~ 8ttM V 2^ 
The complementary error function law describes the tails of the two- 
dimensional eigenvalue support in the region \w\ — 2n ~ O(y'ri). The erfc-law 
can be easily identified as a 'goblet base' in the mean density profile shown 
in Fig. [2] (right panel) . 



• Regime III. — For n 3> 1 and v 1 such that the ratio v/n — q is kept fixed 
(q > 0), we prove that the mean density of complex eigenvalues is described by 
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Figure 4. Destruction of the 'clustering' phenomenon in the regime {n» 1 and 
v 1 with v/n = q kept fixed}. Profiles of the mean density Eq. | |1.32| are plotted 
for n = 100 and q = 1/2 (left panel), q = 1 (middle panel) and q = 3 (right panel). 



the formula 



8tt ^|u>| 2 + n 2 q 2 



. erfc 



\w\ - 2n^fq + l 



(1.34) 



y/2n(q + 2) J ' 

which is an analogue of the Marcenko-Pastur law for complex valued eigenvalues. 
Far away from the critical point \w c \ — 2n^/q + 1, the above equation simplifies to 

1 1 



47T y/\w\ 2 + V 2 



(1.35) 



Notice that (formally performed) v — > limit reduces Eq. (1.34) to Eq. (1.331 as 
expected. Figures [3] and [4] show that the phenomenon of clustering of complex 
eigenvalues around the origin becomes less pronounced with increase of the 
parameter q. 

Having announced the main results of our study, we now turn to their detailed 
derivation. 



2. Non-Hermitean Wishart random matrices at /3 = 2 

2.1. Matrix integral representation of the probability measure 

By definition, the probability measure [§] d7r n-m (VV) induced on \V is provided by a 
two-matrix integral 

*r„, m (W) = J dPi%(X) J dP(%(y) 8 (w - *yt) DW, (2.1) 



where the flat measure DW is defined in Eq. (1.21). The integration in Eq. (2.1) 



can be performed in two different ways leading to two different (albeit equivalent) 
representations of d7r„ !m (VV). 

2.1.1. Integral over complex matrix (Type I) 



We start the evaluation of the two matrix integral Eq. (2.1) by making use of the 
matrix integral representation 

1 



S(A) = 



(2TT) 2n2 



Dq cxp 



tr (qA^ + q^A 



(2.2) 



§ From now on, the notation d7r nim (W) is kept for the probability measure associated with non- 
Hcrmitean Wishart matrices. 
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of the (5-function of a general n x n complex valued matrix A, where the flat measure 
Dq is 



a—l a' — l 



Setting 



a = w xy 



(2.3) 
(2.4) 



we rewrite Eq. (2.1 1 in the form 
<frr„ >m (W) = . \ 2n 5 / Dq exp 



tr gW +q T W 



dP^(Af) / dP^C^exp 
By virtue of the identity 



~tr(^Aft+gt^yt 



DW. 



d^CV) exp 



itr(By + cy 



= exp 



-tr^C) 



(2.5) 



(2.6) 



with the matrices B and C of the size m x n and n x m, respectively, being set to 



B = Cft = ~Xlq, 

(2) 

the integral over dPnJniy) can be performed to bring a reduced representation 



(2.7) 



G?7Tn,m(VV) 



(27T) 



2n 2 



£>q exp 
dP« (AT) exp 



4a 

Applying the (Gaussian integral) identity 



exp 



trlX J DX 



dctr 



DW. 



1„ + -D 

«2 



where the n x n matrix _D is set to 
1 



4a 



7 99 



(2.8) 



(2.9) 



(2.10) 



(2) 

we calculate the integral over dPn,m{X) to arrive at the main result of this Subsection: 



dir n , m {yV) 



1 



(27r) 2 " 2 
x del"™ 



1, 



Dq exp 
1 



tr ( qW + <Z f W 



T qq* DW. 



, ,., , (2.11) 

This matrix integral will be referred to as the type I representation of the probability 
measure d7r„ jm (W) associated with non-Hermitean Wishart random matrices. 

2.1.2. Integral over positive definite Hermitean matrix (Type II) 

To derive an alternative representation of the probability measure d~K n , m (W), we start 
with Eq. (2.8) rewritten in the form 

1 



dir n , m {W) 



(2tt) 2 " : 
x exp 



Dq exp 



tr (gW 1 + q^W 



Vtr (qHX Xt)q 



4a' 



DW. 



(2.12) 
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Now we wish to integrate out the matrix q. In the Wishart domain m > n, this can 
be achieved with the help of yet another easy-to-prove Gaussian integral identity 



D q exp 



i tr (qE 



nt 



Tr™ det~"Fexp 



exp [-tr (gfgt))] 
tr (f^F -1 ^ 



(2.13) 



that holds for arbitrary n x n matrix F and a non-singular, positive-definite n x n 
matrix F. Setting 

1 



2 KV ' 



AT AT 



(2.14) 



we reduce Eq. (2.12) to 



dTT n , m (yV) = 



7T 

x exp 



dP(%(X)det- n (X^) 
-a' 2 tr(w\xX^)- 1 wX\ DVV. 



To facilitate integration over dPn)n{X), we insert an integrated 8- function [jj] 
DQS(Q - XX 1 ) = 1 



Qt=Q 



(2.15) 
(2.16) 



into Eq. (2.151 to come down to 



dn n , m (yv) = (—)' 



i v n 



x exp 



Q f =Q 



L»Q det""Q exp [-a 2 tr Q] 



a' 2 tr (Wq^Vv) / Z?AT <$(Q — XX 1 ) DVV. (2.17) 



Finally, making use of the formula 

[ DX S(Q - XX 1 ) = 7r «(2™-«+i)/2^ j J=L 
•Ao*™ ll 7 =i 



U T=7m 



0(Q) (detQ) r 



(2.18) 



proven in Appendix [Aj we end up with the sought alternative (type II) representation 
of the probability measure d7r„ jm (VV) for non-Hermitean Wishart random matrices: 



d7Tn,m(VV) 



(a 2 a' 2 )"("+™)/ 2 n?="r(i) 



7r n(3n-l)/2 



DQQ(Q)det m - 2n Q exp 



Qt=Q 



Va^tr (Q + VV^Q^VV) 



DW. 



(2.19) 



Equation (2.19) represents the main result of this Subsection. 
2.2. Spectral statistics 

To study the spectral statistics of non-Hermitean Wishart random matrices, we need 
to evaluate the joint probability density function (j.p.d.f.) p n ,m{w\, ■ ■ ■ , w n ) of all 
n complex eigenvalue s of W. Altho ugh both type I and type II matrix integral 
representations [Eqs. (2.11) and (2.19)] of the probability measure d7r„, m (W) can be 
used to this end, the derivation based on Eqs. (2.11 ) is particularly elegant. 

| Here, the flat measure DQ is 



DQ=]JdQ aa JJ dmtQ a:a >£!mQ aiCt/ 



2 Non-Hermitean Wishart random matrices at (3 = 2 



13 



2.2.1. Joint probability density junction of all eigenvalues 



On general grounds, the probability measure dp nm (w\, ■ ■ ■ ,w n ) in the complex eigen- 
value space is related to the probability measure dir n ^ m (W) via the integral transfor- 
mation 



dp ntm {wi, . . . ,w n ) 



[ [ dn n . m (vV = U(w + R)U*) 



(2.20) 

/C/SU(n) JRtC" 

This prescription is a consequence of the Schur decomposition of the n x n complex 
matrix VV with distinct eigenvalues |^] 

VV = U{w + R)U\ (2.21) 

where U is a unitary matrix, R is strictly upper triangular, 

R jk e C, j < k, , . 

o, j > k, (2 - Z2) 

and w = diag(w 1; . . . ,w n ) is a diagonal matrix composed of n complex eigenvalues 
of VV. This decomposition is unique if we label the eigenvalues and require the first 
nonzero element in each column of U to be positive. It holds that 

DVV = DW DRdp,{U), 



(R)jk 



where 



DW = \ A n (w)\ 2 Yldmewj dJm', 



'31 



3 = 1 



DR = 



Y\_ dd\cRjkd3mRjk, 

j<k=l 

and dp(U) denotes the Haar measure on U(n). 



(2.23) 

(2.24) 
(2.25) 



Below, the integrals in Eq. (2.20) will be evaluated based on the type I representa- 
tion 



dTT n , m (W) 



1 



(2tt) 2 ™ ; 
x det" m 



Dq exp 
1 



tr 



q'W 



-qql DW 

4a 2 a 2 



(2.26) 



of the probability measure c?7r njrn (VV). We shall proceed in three steps. 



Step t] 1. — Substitute the Schur decomposed matrix W [Eq. (2.21)] into Eq. (2.26) 
to observe that 

. vol[U(n)l 
dp ntm {wi,...,w n ) = ^tt) 2 " 2 



D£ exp 



tr (|u; t + $}w 



det" 



In + T^-T^ 

4a 2 a 2 



DR exp 



tr ( + ^R 



Here, (, = U^qU is a new integration matrix, and 



vol[U(n)] = / dp(U) 



7r n(«-l)/2 



(2.27) 



(2.28) 



1 See Appendix 33 in Mehta (2004). 
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is the volume of the unitary group U(n). Notice that the integral over 'radial' degrees 
of freedom DR has factored out. 



Step \ 2. — Perform the radial integral 



DR exp 



ReC n( n -l)/2 



tr ( + ^R 



j-<fc=i 



(2.29) 



to find out the remarkable formula 

^n.mK, • • • , «>») = (27r)n( n+l) 



^, e C"<"+ 1 >/ 2 



D'j/j exp 



tr ( i/m^ + i/vu; 



dct~ 



tr, 



4a 2 a 2 



where i/> is an n x n complex valued lower triangular matrix, 
ClA - I °' J < fc > 



(2.30) 
(2.31) 



Step t] 3. — Use Lemma 3 of Appendix |A| to eventually derive: 



2«(a 2 a! i )«(™+ 1 )/ 2 Y\7=i ^) 



T n r>\ 



DW. 



x H\ Wj \ m - n K m _ n (2\ Wj \^W 2 

3 = 1 

Considered together with Eq. ( |2.24 1, this completes our calculation of the j.p.d.f. 

pn,m{wi, ...,«)„) 



(2.32) 



2"(a 2 a 2 )"(™ +i )/ 2 njirr(j) 



T n n\ 



n- =1 r(j)nr=i r (i) 



x|A„(™)| 2 11 \ Wj \ m - n K m . n (2^1 (2.33) 

3=1 

of all n complex eigenvalues for a non-Hermitean Wishart random matrix model defined 
by Eqs. (1.171 - (1.19) in Section 1.2.1 Equation (2.33) represents the main result 
of this Subsection. It can be viewed as a non-Hermitean counterpart of Eq. (1.11) 
considered at f3 = 2. We note that a closely related result was first obtained by 
Osborn (2004) in the context of QCD physics, see also Akemann, Osborn, Splittorff 
and Verbaarschot (2005). 

2.2.2. Correlation functions 

The p-point correlation function is defined in the standard manner (Mehta 2004) 
n! 



R {p) ( 



Wi,- ■ ■ ,w p ) 



(n-p) 



n 

j=p+i 



dSXzwj d3mwj p n ,m(wi, • • • , w p ,w p+1 , ■ ■ ■ ,w n ). 

(2.34) 



Here, 



and 



Pn,m(wi,- ■-,W n ) = C„]„ \A n (w)\ 2 \Wj\ v K v ( 

3 = 1 

n,=i r u) 



(2.35) 



(2.36) 
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The parameter v is defined in Eq. (1.26). [In order to simplify notation, we have set 
a 2 a' 2 = 1/4 in Eg. p33||]. 



To integrate out n — p eigenvalues in Eq. (2.34), we employ the orthogonal polyno- 

(2.37) 

(2.38) 



mials technique (Mehta 2004). Introducing a set of monic polynomials 

Pj(w) = w° + axW 3 ~ + • • ■ 
orthogonal on C with respect to the measure 

dn(w) = \w\ v K v {\w\)d£faw dZmw, 
that is, 



dn(w)pj(w)p k (w) =Nj8j k - 



(2.39) 



Rotational symmetry of the weight function in the measure dfi(w) suggests that Pj (w)'s 
are merely the monomials 

Pj(w) = w 3 
with the normalisation constant Mj being [f] 

Afj = f dfi(w) \w\ 2j = 2 2j+t/+1 7rr(j + l)T(j + v + 1). 
Jc 

This observation allows us to rewrite the j.p.d.f. Eq. ( 2.35[ ) in the form 



(2.40) 



(2.41) 



pn,m{Wx, ■ ■ ■ ,W n ) 



■ dot 



p k -l{Wj) 



det 



Pk~l(Wj) 



fe-1 



nK-r^(K-i)- 



Equivalently, 



where 



Pn,m{w\,'-- >W n ) = ^det[K n , m (lVj,W k )], 

/ ,J\k 



(2.42) 
(2.43) 



k=0 



2 2k k\(k + v)\ 



is the scalar kernel that obeys the projection property 

cKHe £eDm£ K n , m (w, £) K„, m (£, w') = K„, m (w,w')- 



(2.44) 
(2.45) 



By virtue of the Dyson integration theorem (Mehta 1976), Eqs. ( |2.34[ ), ( |2.43[ ) and ( |2.45 l 
imply that the p-point correlation function equals 

Rn.ln(wi,---,W p ) = dctlKn^ilVjiWk)]^^ (2.46) 

as was anticipated. Note that a closely related result appears in the recent work by 
Akemann, Phillips and Shifrin (2009). 



f Notice that 



0=0 



3 = ^nT 



which implies, through the Andreief-dc Bruijn integration formula (Andreief 1883 and de Bruijn 1955) 
that the probability measure Eq. l |2.35| l is properly normalised: 



fl/c 



diKe Wj cDmwj p n ,m(w\, ■ ■ ■ , w„) 



-^-det(/ dftcwd3mww : >w k \w\"K v (\w\) 



n ^ = 1 



3=0 
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2.3. Large-n /large-v analysis of the mean density of eigenvalues 



Equations (|2.44|) and (|2.46|) show that the mean density of complex eigenvalues is 

(2.47) 



determined by the formula 
where the function T n ^(g) equals 



1 



fc=0 



k\{k + v)\ \2 



2k 



(2.48) 



Equation (2.47) is exact as it is valid for arbitrary m > n. Below, wc shall be interested 
in its large-n/large-i' analysis. 



2.3.1. The case n fixed and v ^> 1 (Regime I) 

This region of parameters corresponds to meticulously sampled time series. To study 



the asymptotic behaviour of the mean density of eigenlevels [Eq. (2.47)] in this regime 



we focus on the function J- n ,v{o) defined by Eq. (2.48) to realise that, for v ^> 1, the 



inverse factorial \(k + 1 therein can be approximated by the leading term of the 
Stirling formula, 

1 1 e u 1 



(k + v)\ 

Consequently, we derive 



1 e" ^ 1 



k=0 



2tt v u+1 ' 2 v k " 



tt_ 



v > 1. 



1 e e e ' iv r [ „ 



(2.49) 



(2.50) 

This result prompts that for the large-z^ limit of the spectral density to be well defined, 
the energy variable \w\ must scale with ^/v so that their ratio 



(2.51) 



is kept bounded. Having appreciated this fact, we make use of Lemma 4 of the 
Appendix [B| to treat the modified Bessel function K v {\w\) in Eq. (2.47). As the result, 
Eq. (2.47) boils down to the elegant expression 



> r I n. 



(2.52) 



Equation (2.52) holds for n finite and v 3> 1. A detailed discussion of the qualitative 



behaviour of the mean eigenvalue density Eq. (2.52) has been presented in Section 
1.2.2 see also Fig. [I] (right panel). 



2.3.2. The case n» 1 and v > fixed (Regime II) 

For v ~ O(n ), two different large-n subcases should be distinguished. 



2 Non-Hermitean Wishart random matrices at (3 = 2 



17 



• If |ui| ~ 0{n ), the upper summation limit in Eq. (2.48 ) can be extended to infinity 
yielding g] 

#«(™) = ~/„(HWM)> (2.53) 

where denotes the modihed Bessel function of the first kind. In particular, 

this simple expression is useful in the analysis of the large- n mean density i?i(|u>|) 
both close to the origin (\w\ <C 1) [I] and far away from it (\w\ ^> 1). 



Close to the origin (\w\ <C 1), we obtain 



i , i 



1 

47Ti/ ' 



v = Q, 



(2.54) 



In the opposite limit, the unusual power-law decay of the mean density is observed 

HI 



1 

H»i 47r|ui| 



(2.55) 



Both laws [Eqs. (2.54| and 1 2.55l] are in sharp contrast with the mean density of 
complex eigenvalues in the GinUE (Ginibre's unitary ensemble), a close relative of 
the non-Hermitean Wishart matrix model. In the former ensemble, both domains 
(\w\ <C 1 and \w\ ^> 1) are characterised by the constant mean density (Ginibre 
1965) 

E (l,GmUE), , 1 (2 56) 

7T 

indicative of the emerging Girko law (Girko 1984, Girko 1986, Bai 1997). 



The slow decay at infinity exhibited by Eq. (2.55) implies that this limiting law 

X We have used the Taylor series 



/„(*) 



E 



k\ (k + v)\ \2> 



[See http : //functions . wolfram . com/03 .02.02. 0001 . 01] . 

§ This spectral region is of particular interest in QCD physics, see e.g. Akemann, Phillips and Shifrin 
(2009). 

|| In addition to the Taylor series for I v (z), we have also used the small-z expansion 



;(|)-7, v = o, 



K v (z) 



2 



exhibiting the main terms as z — > 0. Here, 7 : 
http : //functions . wolfram . com/03 . 04 . 06 . 0044 . 01] . 

1 For |arg(z)| < n/2 and \z\ — > 00, it holds: 



and 



I„(z) ■■ 



K„(z) ~ 



12-kz 



0.577215665 . 



l + O 



is the Euler constant. [See 



l + O 



[See http://functions.wolfram.com/03.02.06.0006.01 and . . /03 . 04 . 06 . 0010 . 01.' 
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must break down at some 'critical' point \w\ = \w c \. The mean density normali- 
sation 



2vr 



d\w\ \w\ x 



1 

47r|u;| 



suggests that 



2n. 



(2.57) 



(2.58) 



The position of the break-down point Eq. (2.58) is clearly seen in Fig. [2] (right 
panel). 

Remark. — We note in passing that the algebraically decaying large-n density 
Eq. (2.55) as well as the position of the critical point |w c | = 2n can alterna- 
tively be derived within the two-dimensional Coulomb fluid approach (Chau and 
Zaboronsky 1998, Zabrodin and Wiegmann 2006). Indeed, given the j.p.d.f. of 
the form 



pniwi 



,w n ) oc |A„(u;)| 2 Y[e 



-V(v 



i) 



(2.59) 



it is straightforward to realise that, in the leading order in 1 /n, the mean eigenvalue 
density 



E«H = — V{w,w) 

is proportional to the Laplacian 

d 2 



dwdw 



(2.60) 



(2.61) 



of the confinement potential V{w 1 w). For confinement potentials of the form 
V(w,w) — V(\w\), the eigenvalue density is supported within the disk of the 
radius R — \w c \, where |iu c | is a (positive) solution of the equation 



R 



dV{R) 
dR 



dV{r) 



limr- r 

i — s-o or 



(2.62) 



It is a simple exercise to verify that the large- 1 if | expansion of the modified Bessel 
function if„(|ui|) in Eq. (2.35) yields the effective confinement potential of the 
form 

V eS (\w\ » 1) = \w\ - [v - log M + Odw]- 1 ). (2.63) 

Applying Eqs. (2.60) and (2.62), one reproduces Eqs. (2.55) and (2.58), respec- 
tively. ■ 



• To probe the mean eigenlevel density in the vicini ty of t he cr itical point \w c \ ~ 2n, 
we have to analyse the exact formulae Eqs. (2.47) and (2.48) for \w\ ~ ©(n 1 ). We 
found it useful to put forward the multiplicative ansatz 

R$ m (w) = R&\w)-T c (\w\), (2.64) 

where the function T c (|w|) = T c (|w|;n, v) is to account for (anticipated) signifi- 
cant deviations of the mean density of eigenlevels from the 'naive' limiting curve 
Roo (w) given by Eq. (2.53) in the vicinity of the critical point |w c | ~ 2n. In the 
region 1 <§; \w\ <C \w c \, the function T c (|w|) is expected to approach unity. 
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To determine T c (g) for g £ R + , we first spot that yet another function F nill {g) 
[Eq. ( 2.48| )] entering the exact solution Eq. (2.47), satisfies the differential equation 
2;y 4- 1 1 / n\ 2(n-l) 

KM + K,M ?nAQ) = - v< yri Z S (!) ' t 2 ' 65 ) 

g 1 (n) L (n + is) \ 2/ 

Owing to the relation 

FnAo)= (f\ ^(e)Tc(e), (2.66) 

the function T c (g) is seen to satisfy the differential equation 



l + 2g 



I'M 
h{g) 



K{g) = 



1 



T{n)T{n + v)I v {g) \2 



(2.67) 

By derivation, this equation is exact as it is valid for arbitrary v > 0, n and 
As soon as we are interested in the large-n analysis of a vicinity of the critical 
point \w\ ~ |iu c | ~ 2n, the large-£> reduction of Eq. (2.67) matters. For g >• 1, the 
above equation simplifies to 

T^'(e) + 2 T^(g) - -c n , u g^+v-W e~e, (2.68) 

where 

xFk e 2n 1 

(2.69) 



n ' u 2 2n +»- 5 / 2 r{n)T{n + v) y/nn (2n) 2n +»- 3 / 2 ' 
Focussing on an n Q -vicinity of the critical point by setting 

\w\ = g = 2n + 2n a t, t ~ O(n ), (2.70) 
where the exponent a is yet to be determined, we introduce the new function 

A c (t) = T c (g = 2n + 2n a t; n, v) (2.71) 

which appears to have a well-defined large-n limit provided a = 1/2. This limit is 
described by the equation 



1 



Ai(i) = ~pe 



(2.72) 



Its solution, that approaches unity at t —oo (that is, deep in the bulk to the 
left of the critical point) reads: 

A c (t) = Jerfc(t). (2.73) 



Equivalently, 



Tc(e) 



-erfc 



g — 2n 



(2.74) 



2 ~~~ V 2yjn 

Hence, we conclude that the mean density of eigenlevels in the y/n- vicinity m of 
the critical point is given by the formula 

1 „ {\w\-2n 



8it\w\ 



erfc 



2Jn 



(2.75) 



Notice that Eq. (2.75) predicts that the mean density of eigenlevels at the critical 



point is two times smaller as compared to the naive result Eq. (2.55) if extended 
down to \w\ = |iu c |. This nonperturbative piece of the mean density is responsible 
for forming a 'goblet base' in Fig. [2] (right panel). 



Put differently, Eq. \2.75\ stays valid for 



\w\ - 2n 
2^n 



O(n ) 



kept bounded. 
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The overall behaviour of the mean density is well captured by the single formula 



4^H-^^(N)«v(H)erfc 



| to | — 2n 
2v^ 



(2.76) 



as has been discussed in Section 1.2.2, see also Fig. [2] (right panel). 
2.3.3. The case n» 1 and v/n — q > fixed (Regime III) 

This region of parameters corresponds to the Marcenko-Pastur domain discussed in 
Section [l.l| To study the large-n behaviour of the mean density of complex eigenlevels, 
we follow the idea introduced in Section |2.3.2| Specifically, we decompose the mean 
density as 

Ri%(w) = R£Hw)-f c (\w\), (2.77) 

where 



Rg>{w) = ±I v (\w\)K u (\w\) 



(2.78) 



is the 'naive' limiting curve to be corrected (i n a n onpcrturbative way) by the function 
T c (|to|) = T c (|to|; n, v) proven to satisfy Eq. (2.67). 



Although Eq. (2.781 mirrors Eq. (2.531, it furnishes a limiting law which differs from 
Eq. (2.551; the difference derives from the fact that, in the Regime III, both the ener- 
gies corresponding to the spectrum bulk and the parameter v scale linearly in n. The 
latter necessitates the use of asymptotic expansions for the modified Bessel functions 
due to Olver, Lozier, Boisvert and Clark (2010): 



1 



'2-KV 



and 



I V {VZ) £ 

I' v (yz) = 

K v {yz) 
K' v {vz) ~ - 



VT)(z) 



^ (l + z 2 ) 1 ^^ 



U k (p(z)) 



fc=0 



E 

k=0 



V k {p{z)) 



n e -^(«) ~ k U k (p(z)) 
2^ (llz 2 ) 1 / 4 ^ 1 ' " fe ' 



k=0 



ir e-v^il + z 2 ) 1 ^ 



fe=0 



Here, z£l, and the functions rj(z) and p(z) are defined as 



(z) = \/l + z 2 + log 



n 



p{z) 



i + yr^K 



1 



(2.79) 
(2.80) 

(2.81) 
(2.82) 

(2.83) 
(2.84) 



VT+z 2 " 

Also, U k (p) and V k (p) are polynomials in p of degree 3fc, given by Uo(p) = Vo(p) = 1 
and the recurrence 



dt(l-5t 2 )U k (t), 



(2.85) 



u k+1 {p) = \p 2 ^-p 2 )u' k (p) + \ 

V k+1 {p) =U k+1 (p)- l -p{l-p 2 )U k {p)-p 2 {l-p 2 )U' k {p). (2.86) 



These asymptotic expansions readily bring the large-n formula 

4^ ^J\ w 2\ +V 2 



(2.87) 
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break down at some 'critical' point |«;| = \w c \. The mean density normalisation 

1 



The slow decay at infinity exhibited by Eq. (2.87) implies that this limiting law must 

(2.88) 

(2.89) 



2tt 



d\w\ \w\ x 



47T y/\w\ 2 + V 2 



suggests that 

\w c \ ~ 2ny/q + 1. 
The position of this break-down point is clearly seen in Fig. [4] 



Remark. — Equations (2.87) and (2.89) can alternatively be derived within the two- 
dimensional Coulomb fluid approach discussed in Section [2.3.2| Indeed, making use of 
Eq. (2.81), we observe that the effective confinement potential in the Regime III is 

\/\w\ 2 + v 2 - v log [v + ^\w\ 2 + vA 



K ff (IH) 



M~0(V) 



+-iog(H 2 + ^ 2 ) + o(H- 1 )- 



(2.90) 



Applying Eqs. (2.60) and (2.62), one reproduces, in the leading order in \/v, Eqs. 
(12.87} and (12.89), respectively. ■ 



To describe a vicinity of the critical point \w\ ~ |to c | ~ 2n^Jq + 1, we need to study 
the large- 1 w | reduction of Eq. (2.67). Straightforward calculations based on Eqs. (2.79) 
and ( |2.~80| > yield: 

2v I o 2 - f - 2xV4 

g V v £ 



1 



exp 



-VI) 



(2.91) 

Here, the function r\{z) in the exponent is defined by Eq. (2.83); the constant 
equals 



9(9 + 1) 



^(9+2) 



(2.92) 



^ 2 2 "+^ 5 / 2 r(n)r(n + ^) V 7T ( g + i)n(</+i)(2n)"(9+ 2 )- 3 / 2 ' 
Let us focus on a vicinity of the critical point by setting 

\w\ = q= 2n^/q + 1 + 2n a t, t ~ O(n ), (2.93) 
where the exponent a will be determined later on. Introducing the new function 



hjf) = f c (q = 2n^/q + 1 + 2n a t; n, nqj 



(2.94) 



one observes after somewhat lengthy but straightforward calculations that A c (t) 
appears to have a well-defined large-n limit provided a — 1/2. This limit is described 
by the differential equation 



7T(q + 2) 



exp 



q + 2 



supplemented by the initial condition A c (— oo) = 1. One has: 



Ac(t) = ^erfc ( / 



Equivalcntly, 



T c (g) = -erfc 



q + 2 / 
q - 2n*Jq + 1 

VMq + 2) 



(2.95) 



(2.96) 



(2.97) 
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Consequently, we conclude that the mean density of eigenlevels in the -^/n-vicinity of 
the critical point is given by the formula 



8n Vl< 



| 2 + n 2 q 2 



: erfc 



2n^q + 1 



(2.98) 



Notic e that formally performed q —> limit reproduces Eq. | |2.75 1 derived in Section 
2.3.2|for the case n > 1 and v — O(n ). 



Equations (2.87) and ( 2.98| ) represent the main result of this Subsection. Illustrated 
in Section |1 .2.2 they may naturally be considered as the complex-plane analogues 
of the Marcenko-Pastur law [Eq. (1.13)]. It should be stressed that, contrary to the 
Marcenko-Pastur law for real valued spectra, the mean density of complex eigenvalues 
does not develop a gap around the spectrum origin. 



3. Conclusions 



Driven by potential multidisciplinary applications in statistical analysis of remote com- 
plex systems characterised by distinct sets of stochastic time series, we have introduced 
a non-Hermitean extension of paradigmatic Wishart random matrices and presented a 
detailed study of their spectral properties in the simplest case of complex valued time 
series. In particular, we have determined (i) two alternative one-matrix-integral rep- 
resentations of the probability measure in the space of matrix entries [Eqs. (1.20) and 
(1.22)], (ii) the joint probability density function of all complex eigenvalues [ Eq. (|1.23 l], 
(iii) the eigenvalue correlation functions of arbitrary (finite) order [Eqs. ( 1.24[ ) and 
(1.25)], and also analysed (iv) the mean spectral density in various (large-n/large-^) 
scaling limits [Eqs. ( fT28"| ), ( fl~29| and ( fl~34| )]. 



Interestingly (and also surprisingly for the authors) , a family of closely related random 
matrix models has been introduced, and studied, in a totally different context of the 
QCD physics with a nonvanishing chemical potential (see, e.g., a review by Akemann 
(2007)). This observation (whatever expected or unexpected it may seem) provides one 
more evidence of the mysterious ubiquity of Random Matrix Theory appearing again 
and again in very distant fields of knowledge. 
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Appendices 

A. Three matrix integrals 

In this Appendix, we give a detailed derivation of the three matrix integrals used in 
the main body of the paper. 

Lemma 1. Let Q and S be n x n Hermitean and n x m rectangular complex ma- 
trices, respectively. Then, for m > n, it holds: 

/ DSS(Q- SSA = ^(2m- n +i)/2 ^=i y> ( det Q) m -" 8(Q). (A.l) 

ic-x™ v j n,-=i r u) 



Proof. — Denoting the above integral as I n .m{Q) and employing the matrix integral 
representation 



5 ( A ) = I D(T ex P [* tr ( CT A )l 

of the 5- function of an n x n complex Hermitean matrix A, we set 
A=Q 

to write down: 
1 



In, m (Q) = —-i Dcrexp[itr(crQ)] DS exp f-ttr (sfl(tr - irj)S 

Here and above, the flat measure Dcr equals 



(A.2) 
(A.3) 

(A.4) 
(A.5) 



Dcr = J^J do-jj YJ_ dtyKz o-jk dJmajk, 

j=i j>k=i 

and an infinitesimally small imaginary regulariser —in was introduced to ensure 
convergence of the Gaussian integral 



DS exp 



-itr 



- mj)s)] = {-n:) nm det' m {cr - int n ) . 



Plugging this result back to Eq. (A.4 1, we derive: 

(-i) nm 



On»-n(n-m) / , 



Dcr exp [i tr (crQ) ] det m (cr — ir]t 7 



(A.6) 



(A.7) 



This is the Ingham-Siegel integral of the second type (Ingham 1933, Siegel 1935, 
Fyodorov 2002) calculated in Lemma 2. For m > n, we make use of Eq. (A.9) to 
write down 



/ Dcr exp [i tr (crQ) ] det m (cr — int n ) 

J rr ' rr 



2 „ ngrr(i) 
" u-um 



n(n+l)/2 -nm 



i" m exp[-77trg] (detQ) m -™6(Q). 



(A. 



Substituting Eq. ( |A.8 ) into Eq. (A.7) and letting 77 tend to zero, we reproduce Eq. ( A.l ). 
End of proof. ■ 
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Lemma 2 Q (Ingham- Siegel Integral of the Second Type). Let Q and a be 

n x n Hermitean matrices, and 3m z > 0. Then, for m > n, it holds: 



/ Der exp [i tr (erQ) ] det m (cr - z l r 

J cri —<T 



,nr=Tr(i) 



2" M 3 - 1 t^+wv™ exp ^ tr Q ] (d e tQ) m -™e(g). 



ir=i r c?) 



(A.9) 



Here, Da is given by Eq. (A. 5). 



Proof. — Due to the iterative nature of the forthcoming calculation, it is convenient 
to use the notation Q = Q n and a = a n that highlights the size of matrices in 
Eq. (A.9). Denoting the integral to be treated as F n , m (Q n ), we remark that it may 
only depend on the eigenvalues (Qi, ■ ■ . , Q n ) of Q n . The latter can be written in the 
decomposed form 



Qn i =diag(Q 2 ,...,Q„). 



Qi <-i 

On-l Q n 1 

Accordingly, we represent the matrix er„ in the block form 

Oil 

where an is a real scalar, er„_i is a truncated matrix of the size (n— 1) X (ri — 1), and 



(A. 10) 



(A.ll) 



(Cl2, CT13, ■ ■ • , CTl,n) 



(A.12) 



is a complex vector of the length n— 1. In the above parameterisation, the integration 

(A.13) 



measure in Eq. (A.9) factorises: 



Da n = dan cfiKett n _i cDm«„_i _Dcr„_ 1 
Owing to the identities 

tr (cr„Q n ) = ftnQi + tr (cr n _ 1 Q n _ 1 ), 

and 



(Til - 2 



det (cr„ -zl n ) = det 

V u n-i er?7_l — z Ji„_i 
= det (cr n _! - zl„_i) (cr n - z) - u^^a^i - z t^^ 1 ^^ 
the sought integral F rhm (Q n ) can be decomposed: 



(A.14) 



(A.15) 



CT i-i =cr "-i 



Da n _i exp [i tr (crn-iQ^i) ] det m (cr„„ 1 - z l„_i) 
cDm tt„-i 



x / dcrn exp[icrnQi] 



(en 



J l _ 1 (cr„_i - z\ n -i) 1 u n -i 



(A.16) 



The integral over the scalar an can be calculated by the residue theorem which, for 
3m z > 0, yields 



dan exp[icriiQi] (a n - z) - •u] J _ 1 (cr n _ 1 - zt n -i) 1 u n -i 



= 2tt^—^0(Q 1 )Q™ 1 exp\iQ 1 (z + ul l _ 1 (a n -i-zt n -i) 1 u n -i 



(A.17) 



* Although the proof below closely follows th e one presented by Fyodorov (2002), see also Hua (1963), 
we have included a detailed derivation of Eq. | |A,9[ ) in order to prepare the reader for a proof of Lemma 
3 where essentially the same idea is employed to tackle a more involved matrix integral. 
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Now we are left with the Gaussian integral whose calculation is straightforward: 



(EhXz u n _i 

m— 1 

ITT 

Qi 



cDm u n _i exp 



iQliz + vh^O-n-i-zln-i) 1 «„_i 



det (<r n _i - 2rl n _i) . 



(A.18) 



Plugging this result back into Eq. (A.16), we observe the recurrence relation: 

F n ,m(Q n ) 



n-n+m-lck 



n , Q(Q 1 )QT- n exp[iQ 1 z}. (A.19) 

Assuming that m > n, we iterate the recurrence relation n — 1 times to obtain the 
chain of equations: 

F n — l,m— l \Q n — i) 2 



F n -2,rn-2(Q n -2) 
Fn-2,m-2{Q n -2) 
Fn-3,m-3(Qn-3) 

F 2 (Q 2 ) 
,m— n+1 (Qi) 



r(m - 1) 

2 

T(m- 2) 
2 



7r n - 1 i n+m - 3 e(Q 2 ) Q™"" expfiQaz], 



Jn>— 2 -n+m— 5 



e(g 3 )Qr"ex P [ J Q 3 ^, 



r(m - n + 2) 



2 -m-n+3o 



(A.20) 



where 



Fx^ m - n+ x{Qx) = / dcr„n exp[icr„„Q„] (p nn - z) 
Jr 



m— n+1 



r(m — n + 



r l -m-n+lc, 



e(Q„)Q™-"ex P [iQ„z] 



Combining Eqs. (A.19) - (A. 211 and restoring the initial notation, we derive: 
9 n 

™(n+l)/2.mn exp [j ztr Q] (det Q) m_ "9(Q) . 



F n .m(Q) — -prm p/ ■ \ 

llj'=m-n+l 1 U7 



(A.21) 



(A.22) 



This is equivalent to the r.h.s. of Eq. (A.9|. End of proof. 



(</>) 



Lemma 3. Let ip be an n x n complex valued lower triangular matrix 

0, 3 < k, 
i>jh e C, j > fc, 

and w be an n x n diagonal matrix id = [w\, . . . , w„). It holds that 

i 



(A.23) 



/ 

JiL 



Dt[) exp 



tr I ipw' + tft'w 



dct- m (t n + t/^) 



2n (m -n-l) A" PM IlK-r~ ng m -n(K-|)- 



n™ =1 r(i) ? _, 

Here, K n (w) is the modified Bessel function of the second kind. 



(A.24) 



Proof. — Due to the iterative nature of the forthcoming calculation, it is convenient 
to use the notation if> = ij> n and w = w n that highlights the size of matrices in 
Eq. (A.24). Denoting the integral to be treated as I n ^ n (w n ), we proceed with its cal- 
culation by applying the method used in Lemma 2. Specifically, we decompose the 
n x n matrix xp n as 

' °n-l 



</>„ = 



U n -1 Ipn 



(A.25) 
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where ip nn is a complex scalar, V'n-i is a truncated lower triangular matrix of the size 
(n — 1) x (n — 1), and 



,t 



1 = (^n,li^n,2) ' ' ' , VW-l) 



(A.26) 



is a complex vector of the length n — 1 . In the above parameterisation, the integration 

(A.27) 



measure Dxjj = Dip n in Eq. (A. 24 1 factorises: 



Dip n = eKKe ip nn d3m tj) nn dSKt w n _ieDm m„_i Dip n _ 1 . 
Owing to the identities 

tr (tpnWl + Ip^Wn) = i> nn W n + ilnnWn + tr (V'n-l^l-l + ^n-l^n-l) (A. 28) 

and 



det(l n + ^) = det( W*« T i*U 1+ t^- 1 """ 1 ,/ |2 



det(l„_ x +Vn-l^l-l) 



1 + |VW| 2 + "Li ( !n-l - Vn-l^n-l + V'n-lV'n-l) ^n-l ) 



(A.29) 



the integral I n ,m{ w n) can be decomposed as follows: 



Dip n _ 1 exp 



2 tr (^n-l^l-l + 



x det- m (l„_ 1 + V„_W4_i) / d9M„n / dJm$ 



x exp 



(ipnnWn + 1p nn^n ) 



d3\tu n -i / eDmM„_i 



(A.30) 



(i) Let us hrst calculate the inner integral 



l^nnl) 



1 + |VW| 2 + «Ll ( ln-1 - </>;_l(ln-l + VV-lVd) ) 



t \-X. 



To this end, we introduce a new (vector) integration variable 
where an (n — 1) x (n — 1) Hermitean matrix M n —± is given by 



M n _i 



1 + |Vr, 



(A.31) 
(A.32) 

(A.33) 



Since the flat measure efiKeit rl _i eDmti„_i transforms as follows, 

1 



GKKe u n - 1 cDm « n _i 



detM„_! 

= (1 + |i/'n™|)" _1 det(l n _i +'0 n _ 1 v4_ 1 )c£ftew„_i(0mw„_i ) 



(A.34) 
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the integral Eq. (A. 31) boils down to 

dd\zv n _i I t0m« n _i (1 + «^_ 1 w n _i) (A. 35) 



As the remaining integral is straightforward to calculate 
/ dSKtVn-i l d3mv n -i (l + v] l _ 1 v n -i 



^- i r(T "~" + 1) . (A.36) 
I (m) 



wc conclude that 



L(^„_x, Vl-l, IV^nl) = ^"' ^^^^ (1 + l^l 2 )-" 1 - 1 det(l n _x + tf n _^U). 



(A.37) 



Now, the integral I n ,m{ w n) simplifies to 
,_i r(m — n + 1) 



/n,m(«>Ti) = 



T(m) 



x exp 



■ tr 



d^HeV'nn / GDmi/> n „ exp 



nn^n ) 



(i + Kwn 



2\n-m- 1 



(A.38) 



(ii) The above representation is remarkable as the integral over the reduced matrix 
ip n ~i can be recognised to be / n _i im _i(iu n _i). This observation, considered in 
conjunction with the result 

r i 



dSKttpnn I d3mip nn exp 



yields the recursion relation: 



u/n nn^n ) 



2tt 



2 m -™r(m - n + 1) 



(f + |^ nn | 2 )"- m - 1 
|«g m -™K m _„(K|), (A.39) 



(2tt)» 



In-l,m-l{w n -l) 2™- 1 r(m) 

Assuming that m > n, we iterate it (n — 1) times to obtain the chain of equations: 



(A.40) 



mi— l,m- 


-l(u>n-l) 


(27T)"- 1 




Ml— 2,771 - 


-2(« n _ 2 ) 


2 m - 2 r(m - 


1) 


Ml — 2,771 - 


-2K-2) 


(2tt)"- 2 




*n—3,Tn- 


-3K-3) 


2™- 3 r(m - 


2) 


-*2,m— 71- 


-2(1)2) _ 


(27T) 2 




*l,m— 7i- 


-l(wi) 


2 m -" +1 r(m - 


n 



| Wn _ 1 r-"ir m _„(| U ;„_ 1 |), 



K-2r _n tfm-n(K-2l), 

— li^r-^-nditfal), 



(A.41) 



where 

7i, m - ra+ ] ( «' 1 ! = / r /in ecu / (/JniV'ii exp 

(27T) 1 



(-011^1 + V'nwi) 



(i + i^iir) 



2 \ n — m — 1 



| Wl r-«ir m _„(|t« 1 |). 



2 m -«r(m - n + 1) 
g the 

nKr-"tfm-n(KI). 



Combining Eqs. (A.40) - (A.42) and restoring the initial notation, we derive: 

n n(n+l) 1 
In ' m(W) = 2«( — 1) UT=rn-n +i m f = \ 

This is equivalent to the r.h.s. of Eq. ( A.24[ ). End of proof. 



(A.42) 
(A.43) 
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B. Asymptotic behaviour of the Bessel function K v (x^/v) as v — > oo 

Lemma 4. Let K v (x) be the modified Bessel function of the second kind. The main 
term of the asymptotic expansion of K u (x^/v) as v — > oo equals 

where x > 0. 



Proof. — The easiest way to see that Eq. (B.ll holds true is to start with Basset's 
integral El 



K u (xz) 



1 /2z 

Jn \ x 



r ( v 



i 



dt 



cos(xt) 



( t 2 + z 2)„+l/2 • 



(B.2) 



where 5Hc f > —1/2, a; > 0, and arg(z) < tt/2; the branch of (i 2 + z 2 ) !y+1 / 2 is continuous 
and asymptotic to the principal value of t 2v+1 as < — > +oo. Setting z = y^, we observe 
the large-^ relation 



1 



(t 2 + i/)"+Va 



(B.3) 



Plugging it back to Eq. (B.2 1 and making use of the leading order Stirling formula 



1 



we derive: 



K v {x\/v) 



v/2 



dt cos{xt) e 



(B.4) 



(B.5) 



Performing the integral, we reproduce Eq. (B.ll. End of proof. 
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